Exploring Gas-Phase Chemical Reactions in NH3/B2H6 Systems for Chemical Vapor Deposition Using Reactive Molecular Dynamics.

Hexagonal boron nitride (hBN) has attracted significant attention as a two-dimensional (2D) material due to its unique structure and properties. In this paper, we investigated the gas-phase reactions between B2H6 and NH3 and the reaction pathways potentially leading to the synthesis of hBN by using ReaxFF-based reactive molecular dynamics (MD) simulations. From the hundreds of chemical reaction pathways observed in these MD simulations, we extracted a highly reduced chemical kinetic model to describe the gas-phase mixture evolution in chemical vapor deposition of hBN using B2H6 and NH3 as precursors. The intent is to integrate this chemical model into future computational fluid dynamics (CFD) simulations of actual hBN deposition testing and production reactors to provide enhanced insights for experimental synthesis processes and reactor optimization. The chemical model of this study will serve as a stepping stone for large-scale simulations.

Development and validation of a general-purpose ReaxFF reactive force field for earth material modeling.

ReaxFF reactive force field bridges the gap between nonreactive molecular simulations and quantum mechanical calculations and has been widely applied during the past two decades. However, its application to earth materials, especially those under high T-P conditions relevant to Earth's interior, is still limited due to the lack of available parameters. Here, we present the development and validation of a ReaxFF force field containing several of the most common elements in Earth's crust, i.e., Si/Al/O/H/Na/K. The force field was trained against a large data set obtained from density functional theory (DFT) calculations, including charges, bond/angle distortion curves, equation of states, ion migration energy profiles, and condensation reaction energies. Different coordination environments were considered in the training set. The fitting results showed that the current force field can well reproduce the DFT data (the Pearson correlation coefficient, Rp, is 0.95). We validated the force field on mineral-water interfaces, hydrous melts/supercritical geofluids, and bulk crystals. It was found that the current force field performed excellently in predicting the structural, thermodynamic, and transport properties of various systems (Rp = 0.95). Moreover, possible applications and future development have been discussed. The results obtained in this study suggest that the current force field holds good promise to model a wide range of processes and thus open opportunities to advance the application of ReaxFF in earth material modeling.

Mixing I and Br in Inorganic Perovskites: Atomistic Insights from Reactive Molecular Dynamics Simulations.

All-inorganic halide perovskites have received a great deal of attention as attractive alternatives to overcome the stability issues of hybrid halide perovskites that are commonly associated with organic cations. To find a compromise between the optoelectronic properties of CsPbI3 and CsPbBr3, perovskites with CsPb(BrxI1-x)3 mixed compositions are commonly used. An additional benefit is that without sacrificing the optoelectronic properties for applications such as solar cells or light-emitting diodes, small amounts of Br in CsPbI3 can prevent the inorganic perovskite from degrading to a photo-inactive non-perovskite yellow phase. Despite indications that strain in the perovskite lattice plays a role in the stabilization of the material, a full understanding of such strain is lacking. Here, we develop a reactive force field (ReaxFF) for perovskites starting from our previous work for CsPbI3, and we extend this force field to CsPbBr3 and mixed CsPb(BrxI1-x)3 compounds. This force field is used in large-scale molecular dynamics simulations to study perovskite phase transitions and the internal ion dynamics associated with the phase transitions. We find that an increase of the Br content lowers the temperature at which the perovskite reaches a cubic structure. Specifically, by substituting Br for I, the smaller ionic radius of Br induces a strain in the lattice that changes the internal dynamics of the octahedra. Importantly, this effect propagates through the perovskite lattice ranging up to distances of 2 nm, explaining why small concentrations of Br in CsPb(BrxI1-x)3 (x ≤ 1/4) have a significant impact on the phase stability of mixed halide perovskites.

Watching (De)Intercalation of 2D Metals in Epitaxial Graphene: Insight into the Role of Defects.

Intercalation forms heterostructures, and over 25 elements and compounds are intercalated into graphene, but the mechanism for this process is not well understood. Here, the de-intercalation of 2D Ag and Ga metals sandwiched between bilayer graphene and SiC are followed using photoemission electron microscopy (PEEM) and atomistic-scale reactive molecular dynamics simulations. By PEEM, de-intercalation "windows" (or defects) are observed in both systems, but the processes follow distinctly different dynamics. Reversible de- and re-intercalation of Ag is observed through a circular defect where the intercalation velocity front is 0.5 nm s-1 ± 0.2 nm s.-1 In contrast, the de-intercalation of Ga is irreversible with faster kinetics that are influenced by the non-circular shape of the defect. Molecular dynamics simulations support these pronounced differences and complexities between the two Ag and Ga systems. In the de-intercalating Ga model, Ga atoms first pile up between graphene layers until ultimately moving to the graphene surface. The simulations, supported by density functional theory, indicate that the Ga atoms exhibit larger binding strength to graphene, which agrees with the faster and irreversible diffusion kinetics observed. Thus, both the thermophysical properties of the metal intercalant and its interaction with defective graphene play a key role in intercalation.

Effect of Fe-O ReaxFF on Liquid Iron Oxide Properties Derived from Reactive Molecular Dynamics.

As iron powder nowadays attracts research attention as a carbon-free, circular energy carrier, molecular dynamics (MD) simulations can be used to better understand the mechanisms of liquid iron oxidation at elevated temperatures. However, prudence must be practiced in the selection of a reactive force field. This work investigates the influence of currently available reactive force fields (ReaxFFs) on a number of properties of the liquid iron-oxygen (Fe-O) system derived (or resulting) from MD simulations. Liquid Fe-O systems are considered over a range of oxidation degrees ZO, which represents the molar ratio of O/(O + Fe), with 0 < ZO < 0.6 and at a constant temperature of 2000 K, which is representative of the combustion temperature of micrometric iron particles burning in air. The investigated properties include the minimum energy path, system structure, (im)miscibility, transport properties, and the mass and thermal accommodation coefficients. The properties are compared to experimental values and thermodynamic calculation results if available. Results show that there are significant differences in the properties obtained with MD using the various ReaxFF parameter sets. Based on the available experimental data and equilibrium calculation results, an improved ReaxFF is required to better capture the properties of a liquid Fe-O system.

Correction: Self-limiting stoichiometry in SnSe thin films.

Correction for 'Self-limiting stoichiometry in SnSe thin films' by Jonathan R. Chin et al., Nanoscale, 2023, 15, 9973-9984, https://doi.org/10.1039/D3NR00645J.

Aqueous Stability of Metal-Organic Frameworks Using ReaxFF-Based Metadynamics Simulations.

Aqueous stability is a critical property for the application of metal-organic framework (MOF) materials in humid conditions. The sampling of the free energy surface for a water reaction is challenging due to a lack of a reactive force field. Here, we developed a ReaxFF force field for simulating the reaction of zeolitic imidazole frameworks (ZIFs) with water. We carried out metadynamics simulations based on ReaxFF to study the reaction of water with a few different types of MOFs. We also conducted an experimental water immersion test and characterized the XRD, TG, and gas adsorption properties of the MOFs before and after the immersion test. By considering the energy barrier for a hydrolysis reaction, the simulation results are in good agreement with the experiments. MOFs with open structures and large pores are found to be unstable in metadynamics simulations, where the water molecule can attack or bond with the metallic node relatively easily. In contrast, it is more difficult for water to attack the Zn atom in the ZnN4 tetrahedral structure of ZIFs. We also found that ZIFs with the -NO2 functional groups have higher water stability. Discrepancies between the metadynamics simulation and gas adsorption experiments have been accounted for from the phase/crystallinity change of the structure reflected in the X-ray diffraction and thermogravimetry analysis of the MOF samples.

Step engineering for nucleation and domain orientation control in WSe2 epitaxy on c-plane sapphire.

Epitaxial growth of two-dimensional transition metal dichalcogenides on sapphire has emerged as a promising route to wafer-scale single-crystal films. Steps on the sapphire act as sites for transition metal dichalcogenide nucleation and can impart a preferred domain orientation, resulting in a substantial reduction in mirror twins. Here we demonstrate control of both the nucleation site and unidirectional growth direction of WSe2 on c-plane sapphire by metal-organic chemical vapour deposition. The unidirectional orientation is found to be intimately tied to growth conditions via changes in the sapphire surface chemistry that control the step edge location of WSe2 nucleation, imparting either a 0° or 60° orientation relative to the underlying sapphire lattice. The results provide insight into the role of surface chemistry on transition metal dichalcogenide nucleation and domain alignment and demonstrate the ability to engineer domain orientation over wafer-scale substrates.

Effect of Electrode/Electrolyte Coupling on Birnessite (δ-MnO2) Mechanical Response and Degradation.

Understanding the deformation of energy storage electrodes at a local scale and its correlation to electrochemical performance is crucial for designing effective electrode architectures. In this work, the effect of electrolyte cation and electrode morphology on birnessite (δ-MnO2) deformation during charge storage in aqueous electrolytes was investigated using a mechanical cyclic voltammetry approach via operando atomic force microscopy (AFM) and molecular dynamics (MD) simulation. In both K2SO4 and Li2SO4 electrolytes, the δ-MnO2 host electrode underwent expansion during cation intercalation, but with different potential dependencies. When intercalating Li+, the δ-MnO2 electrode presents a nonlinear correlation between electrode deformation and electrode height, which is morphologically dependent. These results suggest that the stronger cation-birnessite interaction is the reason for higher local stress heterogeneity when cycling in Li2SO4 electrolyte, which might be the origin of the pronounced electrode degradation in this electrolyte.

Self-limiting stoichiometry in SnSe thin films.

Unique functionalities can arise when 2D materials are scaled down near the monolayer limit. However, in 2D materials with strong van der Waals bonds between layers, such as SnSe, maintaining stoichiometry while limiting vertical growth is difficult. Here, we describe how self-limiting stoichiometry can promote the growth of SnSe thin films deposited by molecular beam epitaxy. The Pnma phase of SnSe was stabilized over a broad range of Sn : Se flux ratios from 1 : 1 to 1 : 5. Changing the flux ratio does not affect the film stoichiometry, but influences the predominant crystallographic orientation. ReaxFF molecular dynamics (MD) simulation demonstrates that, while a mixture of Sn/Se stoichiometries forms initially, SnSe stabilizes as the cluster size evolves. The MD results further show that the excess selenium coalesces into Se clusters that weakly interact with the surface of the SnSe particles, leading to the limited stoichiometric change. Raman spectroscopy corroborates this model showing the initial formation of SnSe2 transitioning into SnSe as experimental film growth progresses. Transmission electron microscopy measurements taken on films deposited with growth rates above 0.25 Å s-1 show a thin layer of SnSe2 that disrupts the crystallographic orientation of the SnSe films. Therefore, using the conditions for self-limiting SnSe growth while avoiding the formation of SnSe2 was found to increase the lateral scale of the SnSe layers. Overall, self-limiting stoichiometry provides a promising avenue for maintaining growth of large lateral-scale SnSe for device fabrication.

Joint Theoretical and Experimental Study of Stress Graphitization in Aligned Carbon Nanotube/Carbon Matrix Composites.

Stress graphitization is a unique phenomenon at the carbon nanotube (CNT)-matrix interfaces in CNT/carbon matrix (CNT/C) composites. A lack of fundamental atomistic understanding of its evolution mechanisms and a gap between the theoretical and experimental research have hindered the pursuit of utilizing this phenomenon for producing ultrahigh-performance CNT/C composites. Here, we performed reactive molecular dynamics simulations along with an experimental study to explore stress graphitization mechanisms of a CNT/polyacrylonitrile (PAN)-based carbon matrix composite. Different CNT contents in the composite were considered, while the nanotube alignment was controlled in one direction in the simulations. We observe that the system with a higher CNT content exhibits higher localized stress concentration in the periphery of CNTs, causing alignment of the nitrile groups in the PAN matrix along the CNTs, which subsequently results in preferential dehydrogenation and clustering of carbon rings and eventually graphitization of the PAN matrix when carbonized at 1500 K. These simulation results have been validated by experimentally produced CNT/PAN-based carbon matrix composite films, with transmission electron microscopy images showing the formation of additional graphitic layers converted by the PAN matrix around CNTs, where 82 and 144% improvements of the tensile strength and Young's modulus are achieved, respectively. The presented atomistic details of stress graphitization can provide guidance for further optimizing CNT-matrix interfaces in a more predictive and controllable way for the development of novel CNT/C composites with high performance.

Role of Bilayer Graphene Microstructure on the Nucleation of WSe2 Overlayers.

Over the past few years, graphene grown by chemical vapor deposition (CVD) has gained prominence as a template to grow transition metal dichalcogenide (TMD) overlayers. The resulting two-dimensional (2D) TMD/graphene vertical heterostructures are attractive for optoelectronic and energy applications. However, the effects of the microstructural heterogeneities of graphene grown by CVD on the growth of the TMD overlayers are relatively unknown. Here, we present a detailed investigation of how the stacking order and twist angle of CVD graphene influence the nucleation of WSe2 triangular crystals. Through the combination of experiments and theory, we correlate the presence of interlayer dislocations in bilayer graphene with how WSe2 nucleates, in agreement with the observation of a higher nucleation density of WSe2 on top of Bernal-stacked bilayer graphene versus twisted bilayer graphene. Scanning/transmission electron microscopy (S/TEM) data show that interlayer dislocations are present only in Bernal-stacked bilayer graphene but not in twisted bilayer graphene. Atomistic ReaxFF reactive force field molecular dynamics simulations reveal that strain relaxation promotes the formation of these interlayer dislocations with localized buckling in Bernal-stacked bilayer graphene, whereas the strain becomes distributed in twisted bilayer graphene. Furthermore, these localized buckles in graphene are predicted to serve as thermodynamically favorable sites for binding WSex molecules, leading to the higher nucleation density of WSe2 on Bernal-stacked graphene. Overall, this study explores synthesis-structure correlations in the WSe2/graphene vertical heterostructure system toward the site-selective synthesis of TMDs by controlling the structural attributes of the graphene substrate.

Aqueous phase conversion of CO2 into acetic acid over thermally transformed MIL-88B catalyst.

Sustainable production of acetic acid is a high priority due to its high global manufacturing capacity and numerous applications. Currently, it is predominantly synthesized via carbonylation of methanol, in which both the reactants are fossil-derived. Carbon dioxide transformation into acetic acid is highly desirable to achieve net zero carbon emissions, but significant challenges remain to achieve this efficiently. Herein, we report a heterogeneous catalyst, thermally transformed MIL-88B with Fe0 and Fe3O4 dual active sites, for highly selective acetic acid formation via methanol hydrocarboxylation. ReaxFF molecular simulation, and X-ray characterisation results show a thermally transformed MIL-88B catalyst consisting of highly dispersed Fe0/Fe(II)-oxide nanoparticles in a carbonaceous matrix. This efficient catalyst showed a high acetic acid yield (590.1 mmol/gcat.L) with 81.7% selectivity at 150 °C in the aqueous phase using LiI as a co-catalyst. Here we present a plausible reaction pathway for acetic acid formation reaction via a formic acid intermediate. No significant difference in acetic acid yield and selectivity were noticed during the catalyst recycling study up to five cycles. This work is scalable and industrially relevant for carbon dioxide utilisation to reduce carbon emissions, especially when green methanol and green hydrogen are readily available in future.

Recent Advances in 2D Material Theory, Synthesis, Properties, and Applications.

Two-dimensional (2D) material research is rapidly evolving to broaden the spectrum of emergent 2D systems. Here, we review recent advances in the theory, synthesis, characterization, device, and quantum physics of 2D materials and their heterostructures. First, we shed insight into modeling of defects and intercalants, focusing on their formation pathways and strategic functionalities. We also review machine learning for synthesis and sensing applications of 2D materials. In addition, we highlight important development in the synthesis, processing, and characterization of various 2D materials (e.g., MXnenes, magnetic compounds, epitaxial layers, low-symmetry crystals, etc.) and discuss oxidation and strain gradient engineering in 2D materials. Next, we discuss the optical and phonon properties of 2D materials controlled by material inhomogeneity and give examples of multidimensional imaging and biosensing equipped with machine learning analysis based on 2D platforms. We then provide updates on mix-dimensional heterostructures using 2D building blocks for next-generation logic/memory devices and the quantum anomalous Hall devices of high-quality magnetic topological insulators, followed by advances in small twist-angle homojunctions and their exciting quantum transport. Finally, we provide the perspectives and future work on several topics mentioned in this review.

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics.

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

Toward a Mechanistic Understanding of the Formation of 2D-GaNx in Epitaxial Graphene.

Ultrathin 2D-GaNx can be formed by Ga intercalation into epitaxial graphene (EG) on SiC followed by nitridation in ammonia. Defects in the graphene provide routes for intercalation, but the nature and role of the defects have remained elusive. Here we examine the influence of graphene layer thickness and chemical functionalization on Ga intercalation and 2D-GaNx formation using a combination of experimental and theoretical studies. Thin buffer layer regions of graphene near steps on SiC readily undergo oxygen functionalization when exposed to air or a He/O2 plasma in contrast to thicker regions which are not chemically modified. Oxygen functionalization is found to inhibit Ga intercalation leading to accumulation of Ga droplets on the surface. In contrast, Ga readily intercalates between EG and SiC in the thicker graphene regions that do not contain oxygen. When NH3 annealing is carried out immediately after Ga exposure, 2D-GaNx formation is observed only in the oxygen-functionalized regions, and Ga intercalated under thicker nonfunctionalized graphene does not convert to GaNx. Density functional theory calculations demonstrate that oxygen functionalization of graphene alters the binding energy of Ga and NH3 species to the graphene surface. The presence of hydroxyl groups on graphene inhibits binding of Ga to the surface; however, it enhances the chemical reactivity of the graphene surface to NH3 which, in turn, enhances Ga binding and facilitates the formation of 2D-GaNx. By modifying the EG process to produce oxygen-functionalized buffer layer graphene, uniformly intercalated 2D-GaNx is obtained across the entire substrate surface.

Development of a Mg/O ReaxFF Potential to describe the Passivation Processes in Magnesium-Ion Batteries.

One of the key challenges preventing the breakthrough of magnesium-ion batteries (MIB) is the formation of a passivating boundary layer at the Mg anode. To describe the initial steps of Mg anode degradation by O2 impurities, a Mg/O ReaxFF (force field for reactive systems) parameter set was developed capable of accurately modeling the bulk, surface, adsorption, and diffusion properties of metallic Mg and the salt MgO. It is shown that O2 immediately dissociates upon first contact with the Mg anode (modeled as Mg(0001), Mg(10 1 ‾ $\bar 1$ 0)A, and Mg(10 1 ‾ $\bar 1$ 1)), heating the surface to several 1000 K. The high temperature assists the further oxidation and forms a rock salt interphase intersected by several grain boundaries. Among the Mg surface terminations, Mg(10 1 ‾ $\bar 1$ 0)A is the most reactive, forming an MgO layer with a thickness of up to 25 Å. The trained force field can be used to model the ongoing reactions in Mg-air batteries but also to study the oxidation of magnesium metal in general.

Role of tilt grain boundaries on the structural integrity of WSe2 monolayers.

Transition metal dichalcogenides (TMDCs) are potential materials for future optoelectronic devices. Grain boundaries (GBs) can significantly influence the optoelectronic properties of TMDC materials. Here, we have investigated the mechanical characteristics of tungsten diselenide (WSe2) monolayers and failure process with symmetric tilt GBs using ReaxFF molecular dynamics simulations. In particular, the effects of topological defects, loading rates, and temperatures are investigated. We considered nine different grain boundary structures of monolayer WSe2, of which six are armchair (AC) tilt structures, and the remaining three are zigzag (ZZ) tilt structures. Our results indicate that both tensile strength and fracture strain of WSe2 with symmetric tilt GBs decrease as the temperature increases. We revealed an interfacial phase transition for high-angle GBs reduces the elastic strain energy within the interface at finite temperatures. Furthermore, brittle cracking is the dominant failure mode in the WSe2 monolayer with tilted GBs. WSe2 GB structures showed more strain rate sensitivity at high temperatures than at low temperatures.

Importance of Nuclear Quantum Effects on Aqueous Electrolyte Transport under Confinement in Ti3C2 MXenes.

Protons display a high chemical activity and strongly affect the charge storage capability in confined interlayer spaces of two-dimensional (2D) materials. As such, an accurate representation of proton dynamics under confinement is important for understanding and predicting charge storage dynamics in these materials. While often ignored in atomistic-scale simulations, nuclear quantum effects (NQEs), e.g., tunneling, can be significant under confinement even at room temperature. Using the thermostatted ring polymer molecular dynamics implementation of path integral molecular dynamics (PIMD) in conjunction with the ReaxFF force field, density functional tight binding (DFTB), and NequIP neural network potential simulations, we investigate the role of NQEs on proton and water transport in bulk water and aqueous electrolytes under confinement in Ti3C2 MXenes. Although overall NQEs are relatively small, especially in bulk, we find that they can alter both quantitative values and qualitative trends on both proton transport and water self-diffusion under confinement relative to classical MD predictions. Therefore, our results suggest the need for NQEs to be considered to simulate aqueous systems under confinement for both qualitative and quantitative accuracy.

Evaluation of the Reax Force-Field for Studying the Collision of an Energetic Proton with the DNA.

The early DNA damage induced by ionizing radiation depends on how ionizing particles transfer energy to this molecule and the surrounding medium, mostly water. In preliminary studies, we found that the energy transferred by a 4 keV proton to a cytosine-guanine base pair in a classical simulation collision using the ReaxFF potential is much smaller than that obtained by a quantum calculation using time-dependent density functional theory (TDDFT). We observed that there are two main reasons for that: no accurate force-field for this situation and problems while dealing with the proton charge during the collision. Here, we only focus on the interaction potential. We calibrated the van der Waals energy term of the ReaxFF potential using TDDFT calculations and a genetic algorithm, specifically for the interaction of a proton with the DNA constituent atoms (carbon, hydrogen, phosphorus, nitrogen, and oxygen). We obtained a significant improvement in the interaction potential and, consequently, in the scattering angle of the proton colliding with the target atoms in question. However, we conclude that despite the improvement for the force-field and scattering angle, the classical charge equilibration method should also be improved to properly describe the proton-DNA collision process.